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At a distance of 249 Mpc (z=0.056), Cygnus A is the only powerful FRII radio galaxy for which a 
detailed sub-parsec scale imaging of the base of both jet and counter-jet can be obtained. Observ¬ 
ing with VLBI at millimeter wavelengths is fundamental for this object, as it uncovers those 
regions which appear self-absorbed or free-free absorbed by a circumnuclear torus at longer 
wavelengths. We performed 7 mm Global VLBI observations, achieving ultra-high resolution 
imaging on scales down to 90 /ias. This resolution corresponds to a linear scale of only ^400 
Schwarzschild radii. We studied the transverse structure of the jets through a pixel-based analysis, 
and kinematic properties of the main emission features by modeling the interferometric visibilities 
with two-dimensional Gaussian components. Both jets appear limb-brightened, and their opening 
angles are relatively large (0j ~ 10°). The flow is observed to accelerate within the inner-jet up to 
scales of ~ 1 pc, while lower speeds and uniform motions are measured further downstream. A 
single component seen in the counter-jet appears to be stationary. These observational properties 
are explained assuming the existence of transverse gradient of the bulk Lorentz factor across the 
jet, consisting of a fast central spine surrounded by a slower boundary layer. 


12th European VLBI Network Symposium and Users Meeting, 
7-10 October 2014 
Cagliari, Italy 


* Speaker. 

^Member of the International Max Planck Research School for Astronomy and Astrophysics at the Universities of 
Bonn and Cologne. 


© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 


http://pos.sissa.it/ 








High resolution mm-VLBI imaging of Cygnus A 


Bia Boccardi 


1. Introduction 

The powerful radio galaxy Cygnus A (z=0.056, Dl= 249 Mpc) is a unique target for investi¬ 
gating the physical processes taking place at the very base of a relativistic jet. According to current 
magneto-hydrodynamic models and simulations [1], the crucial mechanisms for jet launching and 
acceleration happen on very small scales, of the order of tens to hundreds Schwarzschild radii, 
which can be reached by present-day mm VLBI arrays when observing the nearest sources. Global 
VLBI observations^ of Cygnus A at 7 mm (43 GHz) enable high resolution imaging on scales down 
to ^90 jUas, translating into a linear size of only 400 Schwarzschild radii at the source redshift (for 
4 ^bh = 2.5 X lO^M© [2]). The combination of the small observing beam and the high radio fre¬ 
quency provides a sharp view towards the central engine and unveils those regions which appear 
either self-absorbed or that are affected by free-free absorption from a molecular torus at lower 
frequencies (e.g. [3]). Showing a prominent counter-jet also on parsec scale, Cygnus A is likewise 
ideal for testing the unified scheme of AGN [4], with the additional advantage that both the jet and 
the counter-jet can be transversally resolved. While spatial resolution is essential for shedding light 
on the aforementioned topics, such detailed imaging pays back by revealing a very complex picture 
of the internal physics occurring in the plasma flow. In particular, there is increasing evidence that 
jets are not homogeneous outflows, but can show a pronounced transverse stratification. For exam¬ 
ple, limb brightened jets have been observed on pc-scale both in some nearby FR I radiogalaxies, 
like M 87 [5] and 3C 84 [6] and in blazars, as in the case of Mrk 501 [7] and 3C 273 [8]. This feature 
has been interpreted with the existence of a spine-sheath structure in the jet, comprising a central, 
light and ultra-relativistic spine surrounded by a layer of slower and denser material [9]. Whether 
this structure arises directly from the jet formation process or is instead the result of instabilities in 
the flow, is still unclear. Anyhow, it is important to keep in mind that, because of relativistic effects, 
the presence of a transverse gradient of the bulk Lorentz factor can strongly affect the observational 
properties of jets. This is especially true for sources like Cygnus A whose jets are seen at a large 
viewing angle. In this case, the fast part of the flow can be completely invisible, both in jet and 
counter-jet, because of relativistic de-boosting (Fig. 1). In the following we present a study of 
the jet transverse structure and kinematics in Cygnus A from Global VLBI observations at 7 mm, 
showing that Doppler de-boosting may indeed play a major role. 



Figure 1: Doppler factor d versus (large) viewing angle 6 for different intrinsic speeds /3. 


^The European VLBI Network is a joint facility of European, Chinese, South African and other radio astronomy institutes funded 
by their national research councils. The National Radio Astronomy Observatory is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 
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Table 1: Log of observations and map characteristics. Col. 1: Day of observation. Col. 2: VLBA - Very Long Baseline 
Array; GBT - Green Bank Telescope; On - Onsala; Nt - Noto; Eb - Effelsberg; Yb - Yebes. Col. 3: Beam and position 
angle. Col. 4: Peak flux density. Col. 5: Map noise. All values are for untapered data with uniform weighting. 


Date 

Antennas 

Beam EWHM [mas, deg] 

Speak [mJy/beam] 

rms[mJy/beam] 

23/10/2007 

VLBA, GBT, Yb, Eb, On, Nt 

0.23x0.11,-21.9° 

176 

0.13 

16/10/2008 

VLBA, GBT, Yb, Eb, On, Nt 

0.22x0.10,-11.0° 

256 

0.11 

19/03/2009 

VLBA, GBT, Yb, Eb, On, Nt 

0.23x0.09,-12.5° 

264 

0.14 

11/11/2009 

VLBA, GBT, Yb, Eb, On, Nt 

0.19x0.10,-20.1° 

235 

0.32 


2. Data set 

Our Cygnus A data set comprises four epochs from Global VLBI observations made between 
2007 and 2009 (see Table 1). For the first three epochs, data were recorded in dual polariza¬ 
tion mode using 4 subbands (IFs) with a total bandwidth of 64 MHz per polarization (512 Mb/s 
recording rate). The last epoch was observed in single polarization, with 16 subbands and a to¬ 
tal bandwidth of 128 MHz. The data were calibrated in AIPS following the standard procedures, 
while imaging and self-calibration of amplitude and phase were performed in DIFMAP. In Fig. 2, 
we present an image obtained by stacking the CLEAN images from the four epochs, after restor¬ 
ing them with a common circular beam of 0.1 mas. The source shows at this frequency a complex 
structure, with an extended counter-jet and a prominent limb-brightening of the jet on mas scales. 



0.00091 0.0056 0.013 0.024 0.039 0.056 0.076 0.1 0.13 

[Jy] 

Figure 2: Stacked 43 GHz image from observations in 2007-2009, rotated by 16° clockwise. It is convolved with a 
circular beam of 0.1 mas EWHM. Contours represent isophotes at 0.3, 0.6, 1.2, 2.4, 4.8, 9.6, 19.2, 38.4, 76.8, 153.6 
mJy/beam. The four maps were aligned as described in Pig. 5. 


3. Ridge line 

We carried out a detailed analysis of the jets transverse structure and opening angle for a 
single epoch (November 2009). The image was restored with a circular beam of 0.15 mas FWHM, 
rotated by 16° and sliced pixel-by-pixel (1 px=0.03 mas) in the direction perpendicular to the jet 
axis. Then, a Gaussian fit of the brightness distribution along the y axis was performed twice in 
each slice, first fitting only one Gaussian and then fitting two Gaussian profiles. We confirm earlier 
findings [10] of a double ridge line on the jet side, and we also find evidence for it in the counter-jet 
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Figure 3: Double ridge line in jet and counter-jet in the 43 GHz map from November 2009. Vertical lines are errorbars 
on the peak positions of the fitted Gaussians, and are set to 1/5 of the FWHM of the respective component. Note that the 
X and j scales are different. 

(Fig. 3). Here the limb brightening is less pronounced than in the jet side, the two rigde lines are 
closer, but still clearly separated after accounting for positional errors, estimated to be 1/5 of the 
FWHM of the respective component. In the outer-jet the northern rail appears to fade, while the 
southern rail is visible up to a core separation —3.5 mas. 

4. Opening angle 

The jet width was calculated for each slice as the separation between the edges of the two 
(deconvolved) FWHM. The results are shown in Fig. 4 (left), where the deconvolved widths ob¬ 
tained from a single-Gaussian slice fit are also plotted for comparison, and show good agreement. 
Both the jet and the counter-jet appear conical, with some minor oscillations of the width, up to 
r ~ zb 1.7 mas, where a prominent narrowing takes place. For \r\ <1.7 mas, the half opening angle 
0^PP was calculated by performing a linear fit to the jet width (Fig. 4, right), which shows for the 
counter-jet a factor of two narrower opening angle than for the jet: 

<^'*PP = (5.10±0.16)°, = (2.42 ±0.21)° (4.1) 

In the classical scenario of intrinsic symmetry of the jet and the counter-jet, their apparent opening 






Figure 4: Left - Jet width versus core separation as deduced from fitting two Gaussian functions (magenta) or a single 
Gaussian (green). Right - Half apparent opening angles calculated from a linear fit to the jet width. 
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angle should, in principle, be the same, unless they are strongly misaligned. A small misalignment 
is indeed present in Cygnus A both on pc and kpc scales, but it is not enough to justify the observed 
difference. An alternative explanation may invoke the difference between the signal-to-noise ratios 
(SNR) of approaching and receding sides: the considerably dimmer counter-jet may indeed appear 
to be narrower than the brighter jet. However, the hypothesis of intrinsic asymmetry cannot be 
readily discarded with the present data. As it was also found in other radio galaxies, the opening 
angles in Cygnus A are considerably larger than what it is found in blazars, typically featuring full 
intrinsic opening angles of (1—2)° [11]. Assuming a reasonably large viewing angle of 70°, the 
intrinsic full opening angle of Cygnus A is ^ 10° on the jet side. This result supports the spine- 
sheath scenario, in which sources seen at larger viewing angles become more and more sheath- 
dominated, and therefore broader, while the emission from blazars mostly comes from the strongly 
boosted, thin spine. 

5. Kinematic structure 

Finally, we present a study of the kinematic properties of the radio emission in the jets during 
the 2-year monitoring interval. In order to obtain a simplified model of the source, describing only 
the main emission features, we have performed a MODELFIT analysis in DIFMAP, i.e we fit 
circular Gaussian patterns to the self-calibrated visibilities (Fig. 5, left). The images were aligned 
to the position of the component N, assumed to be stationary. The same choice was made in [3], also 
based on compactness and spectral index arguments. According to our cross-identification, a new 
component (Jl) was ejected between October 2007 and March 2009, giving rise to a brightening 
of about +50% of the core and inner-jet region. The proper motion of components J2 and J3 
is best described by quadratic functions, indicating an acceleration, while features in the outer- 
jet move each at a roughly constant speed. The apparent speed jSapp = Vapp/c can be calculated 
as jSapp = +z)), where the luminosity distance, c is the speed of light and z is 

the redshift. The dependence of jSapp on the core separation is shown in Fig. 5 (top right). In 
the approaching jet we observe an acceleration up to jSapp = 1.18 ±0.12 within ^ 0.8 mas from 
the core. Then the speed drops significantly, remaining largely sub-luminal in the outer-jet. At 
a core separation r ^ —2 mas, we observe a stationary feature that approximately coincides with 
the position of the narrowing in Fig. 4 (left). Component Cl in the counter-jet does not show 
a significant proper motion (jSapp = 0.004 ±0.060). These characteristics can be well explained 
by the stratification scenario. By looking at the light-curves of the components in Fig. 5 (bottom 
right), we notice that those crossing the acceleration region (Jl, J2, J3) also show a decreasing flux 
density with time, while flux densities are more constant for components in the outer-jet. This could 
be another piece of evidence for the existence of a fast spine, which is accelerating and therefore 
getting fainter, until the emission gets fully dominated by the slower layers of the sheath. This 
also naturally explains the observed limb-brightening. A similar argument can be applied to the 
counter-jet: only the very slow, and therefore the less de-boosted part of the flow can be seen. 

6. Conclusions 

We have presented preliminary results from a high resolution millimeter VLBI study of the 
radio galaxy Cygnus A. We have shown that the observed limb brightening and the kinematic 
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Figure 5: Left - 43 GHz modelfit maps of Cygnus A. Images were tapered and convolved with a common beam 
FWHM of 0.2 X 0.1 mas in P.A. -15°. Contours represent isophote curves at 1.8, 3.6, 7.2, 14.4, 28.8, 57.6 and 115.2 
mJy/beam. The reference point for the alignment is component N, assumed to be stationary. Top right - Apparent speed 
versus core separation. Bottom right - Light curves of the modelfit components. 


properties of the jet, as well as its large opening angles, can be well explained in a scenario of 
transverse jet stratification. The impact of stratification on the observational properties of jets 
depends on the viewing angle, and it is certainly strong for radiogalaxies. 
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